Melting of Discrete Vortices via Quantum Fluctuations 
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We consider nonlinear boson states with a nontrivial phase structure in the three-site Bose- 
Hubbard ring, quantum discrete vortices (or q-vortices), and study their "melting" under the action 
of quantum fluctuations. We calculate the spatial correlations in the ground states to show the 
superfluid-insulator crossover and analyze the fidelity between the exact and variational ground 
states to explore the validity of the classical analysis. We examine the phase coherence and the effect 
of quantum fluctuations on q-vortices and reveal that the breakdown of these coherent structures 
through quantum fluctuations accompanies the superfluid-insulator crossover. 
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Vortices are fundamental objects in physics which ap- 
pear in different fields including phase singularities in 
optics Q and circulating bosons in Bose-Einstein con- 
densates (BECs) 0]. Periodic lattices such as periodic 
photonic structures for light waves or optical lattices 
for BECs can modify strongly the wave propagation, 
and may support novel types of vortex states termed 
discrete vortex solitons 0, U Such discrete vortices 
describe spatially localized circular energy flows. They 
have been studied theoretically and observed experimen- 
tally in optically-induced photonic lattices as stable self- 
trapped states of light carrying a nontrivial angular mo- 
mentum Similar localized vortex states have been 
predicted to occur for BECs in optical lattices 0, . 

Since many of classical lattice systems supporting dis- 
crete vortices originate from quantum models or have a 
well-defined quantum limit, the fundamental question is: 
Do these coherent structures survive under the action 
of quantum fluctuations? In this Letter we answer this 
question and reveal the intimate connection between the 
classical discrete vortex (CDV) and its quantum counter- 
part, which we call quantum discrete vortex or q-vortex) 
In some sense, these q-vortex states can be compared 
with quantum breathers 0, 0, Q which provide a quan- 
tum analog of the self-trapped states with localized energy 
in discrete classical lattices. 

In this Letter, we consider the simplest case of a three- 
site Bose-Hubbard ring. First, reducing the model to 
the discrete self-trapping equations with a variational 
approach, we find CDVs. Next, calculating the spatial 
correlations in the ground states and the fidelity be- 
tween the exact and variational ground states, we study 
the superfluid-insulator crossover and find a valid regime 
of the classical variational approach. Last, analyzing 
the phase coherence and quantum fluctuations of the q- 
vortices, we find that the coherent structure "melts" un- 
der the action of quantum fluctuations, and this melting 
process accompanies the superfluid-insulator crossover. 
This is in a sharp contrast to CDVs, whose phase co- 



herence is independent of the inter-site coupling strength 
and the total number of particles in the vortex. 




FIG. I: (colour online) Schematic diagrams, (a) Bosons in 
a strongly trimerized Kagome lattice |10| can be reduced to 
three-site Bose-Hubbard rings, (b) Combining a proper two- 
dimensional harmonic potential with the Kagome lattice, a 
triple-well potential can be generated; (c) three-site state. 

The basic concepts of the quantum discrete vortices 
can be captured by the simplest case of a three-site Bose- 
Hubbard ring, as shown schematically in the diagrams of 
Fig. 1. The triple- well potential required for such a state 
may be readily found as a subset of the familiar Kagome 
lattices |l(J [see Fig. l(a,b)]. For a deep potential, the 
three-site subsets become decoupled from one another 
and each of them obeys the Hamiltonian, 
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H=-T J2 (a+a m + a l a+) + -Y / n l (n l -l), (1) 

<l,m> 1=1 

where the operators and a; generate and annihilate 
a particle on the Z-site, and T and U stand for the tun- 
neling and on-site interaction strengths, respectively. We 
consider a repulsive inter-particle interaction (U > 0) as- 
suming the application of our results to BECs. 

Our first goal is to connect the classical and quantum 
pictures for describing the stationary states of this sys- 
tem, and to this end we use the time-dependent varia- 
tional approach 0|- F° r the three-site model, we intro- 
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instead of a conventional product of Glauber's coherent 
states. Here, N is the total number of bosons, \vs) is 
the vacuum state \n\ = 0, rii = 0, 7^3 = 0), and the com- 
plex amplitudes ipi satisfy the normalization condition, 
iV'il 2 + I "02 1 2 + 1 1P3 1 2 = iV. In contrast to the variational 
approach based on Glauber's coherent states this 
variational approach conserves the total particle number 
N. Minimizing the corresponding action, we derive the 
classical Hamiltonian for the complex amplitudes ipi, 



H c = -T J2 + + 
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where U N = (N - l)U/N. The extra factor (N - l)/N 
does not appear in the Hamiltonian obtained from the 
variational approach using Glauber's coherent states, and 
the complex amplitudes ipi obey the three-site discrete 
self-trapping equations [!J, [F|, [lj] : 



% dt 



+ u N \^\%, 
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with {l,m} = 1,2,3. We are interested in the sta- 
tionary solutions and therefore make the substitution 
ipi = A\ exp[— i(p± — (pi)] with chemical potential /1, non- 
negative real amplitudes Ai and phases (pi. The sys- 
tem of equations (4) has a well-known degenerate set 
of solutions with equal amplitudes ^J. The solutions 
with cpx = <j>2 = 03, A\ = A 2 = A 3 = y/N/3 and 
/i = (N — l)U/3 — 2T describe the ground state of the 
system. The symmetric vortex states have A\ — A2 = 
A3 = y/N/3, (p2 — <pi — 4>3 — <p2 = 2ln/3 (where I are 
nonzero integers) and /i = (N — l)U/3 — 2Tcos(2/7r/3). 
These three-site CDVs exist for arbitrary values of N, 
U and T. There exist also asymmetric solutions which 
have only two equal amplitudes. These states are the 
excited states of the system and ultimately connect with 
localization maintained by self-trapping. 

We now return to the fully quantum description for 
the Hamiltonian in which the behavior depends on 
the ratio U/T and N. The most prominent effect in in- 
finite quantum systems is the superfluid-insulator tran- 
sition. In our finite-sized system, the ground state is in- 
deed superfluid in the strong tunneling limit, U/T <C 1, 
as expected from the infinite case. In the weak-tunneling 
limit, U/T ^> 1, the ground state is insulating for com- 
mensurate cases (N — 3k with a positive integer k) but 
has a small superfluid fraction accompanying an insulat- 
ing core for incommensurate cases (N ^ 3fc). However, 
instead of a sharp phase transition observed for infinite 
systems, these two limits are connected by a crossover 



regime. That is to say, the superfluid-insulator crossover 
in this finite-site system with commensurate fillings takes 
the place of the superfluid-insulator transition predicted 
and observed in the infinite systems |lfil | , and hence there 
is no well-defined critical value of the ratio U/T between 
the superfluid and insulating phases. Using exact numer- 
ical diagonalization, we calculate the spatial correlations 
Ka£" cij_|_i) J /n aV) where n av = 7V/3 is the average number 
of particles per site, see Fig. 2. These results show that 
in the limit of weak tunneling the superfluid fraction de- 
pends strongly on TV and the spatial correlations decrease 
with U/T such that the system approaches its classical 
counterpart analyzed above. 
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FIG. 2: (colour online) Spatial correlations vs. U/T. 



As the first step to compare the classical results found 
using the variational approach with those obtained using 
exact diagonalization, we compare the ground-state en- 
ergies. Denoting E^ a the ground-state energy obtained 
from the variational approach and E^ x that obtained by 
exact diagonalization, we show that the fractional dif- 
ference of the ground state energies (1 — E^ a /E^ x ) in- 
creases with the ratio U/T close to the linear limit, 
U = 0. In the limit of strong tunneling, i.e. for 
U /T <C 1 , the results show an excellent agreement be- 
tween the two approaches. This is also confirmed by the 
fidelity between the exact and variational ground states, 
F = | va (GS|GS) CX | 2 , which decreases with U/T for arbi- 
trary AT, see Fig. 3. With U/T < 0.1995, for N from 1 
to 18, (1 - Eg/ Eg) < 8 x 1CT 3 and F > 0.9790. The 
fidelity decrease with U/T indicates that the breakdown 
of the variational approximation accompanies the occur- 
rence of the superfluid-insulator crossover. 

Now we can construct the quantum vortex states. The 
fully quantum counterparts for symmetric CDVs can be 
presented as SU(3) coherent states |^) with complex am- 
plitudes Vfc = \/N/3 x exp [i (ip + 2(k - 1)£tt/3)] with 
k(= 1,2,3) and L denote the site index and the vortex 
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FIG. 3: (colour online) Fidelities between the variational and 
exact ground states vs. U/T for different N. 



state. Adiabatically varying U/T to the strongly nonlin- 
ear limit, these q- vortices evolve into different final states 
dependent on the charges L. For L = 3k, the final states 
are single-peaked states in the distribution of the prob- 
ability amplitudes |C(ni, n 2 )| 2 (the ground state). How- 
ever, for L = 3fc ± 1, the quantum vortex (5) appears 
to be a well-defined eigenstate with a desired probabil- 
ity P es if U/T is less than a certain value. For N = 6, 
P es > 0.9850 when U/T < 0.1259. As the nonlinearity is 
adiabatically increasing, the q- vortex breaks down ending 
up in the limit of strong nonlinearity as a triple-peaked 
state in the probability distribution. The appearance of 
the single-peaked and triple-peaked states is independent 
of total atom number, it indicates the loss of the circular 
current of particles and effective melting of the discrete 
vortex structure. In Fig. 4, we show the quantum adia- 
batic evolution of the q- vortices with TV = 6. 



charge, respectively. Such q-vortices, 

N N-n 2 

|QDV) = Y, Yl C («2,n 3 )|n 1 ,n 2 ,n 3 ), (5) 

712=0 ^3 — 

can be prepared in the limit of strong tunneling 
by imprinting the specific phase structure (</?, ip + 
2Lir/3,tp + 4Lir/3) onto the ground states, |GS) = 
53 X) Gr(ri2, TI3) \ni,Ti2, TI3), and the coefficients satisfy 

C(n 2 , n 3 ) = G(n 2 ,n 3 ) exp U[N(p + ^n 2 Lir + ^n 3 Lir]^ , 

where N = ni + n 2 + n 3 . Fast phase imprinting for an 
atomic cloud with a required phase structure via atom- 
laser interaction was also suggested for entanglement 
preparation as wei l as atomic dark-soliton gener- 
ation and measurement [lij . The asymmetric q-vortices 
can be obtained by applying similar phase imprinting 
procedures to excited states whose classical limit (strong 
tunneling limit) correspond to asymmetric solutions. Be- 
low, we consider only symmetric vortices. 

In the linear limit, the q-vortices with L = 3fc (fc are 
integers) are ground states of the system, while those 
with L = 3fc ± 1 correspond to the excited states with 
higher energies. Beyond the linear limit, the quantum 
counterparts for CDVs can be obtained by the quantum 
adiabatic evolution. For our small-size systems we simu- 
late the time evolution in the complete Hilbert space with 
the Runge-Kutta integration scheme pjj . We change the 
parameters adiabatically such that the populated state is 
always very close to an eigenstate. To determine whether 
a given state or even the quantum vortex itself is a well- 
defined eigenstate, we project it onto a complete Hilbert 
basis of an ensemble of the orthogonal eigenstates. In the 
linear limit, we find that all nonzero-probability compo- 
nents of the vortex have identical eigenvalues and, there- 
fore, we may conclude that the q- vortex is an eigen- 
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FIG. 4: (colour online) Quantum adiabatic evolution of q- 
vortices with TV = 6 from linear to nonlinear limit, (a) and (b) 
correspond to the cases of L = 3fc and L = 3fc±I, respectively. 

To explore in detail how the classical discrete vortices 
"melt" , we calculate quantum fluctuations of their quan- 
tum counterparts, i.e., the q-vortices. Using the quan- 
tum phase concept | 20| , we introduce the cosine and sine 
functions of the quantum phase 4>j for the j—tYi site as 
cos tpj = ^ (a J + Oj ) and sin <f>j = (a J — a,j ) with the 
constant Kj determined by the particle number. Thus, 
the cosine and sine functions for the two-body phase dif- 
ference <j) — (f>2 — 4>i can be defined as, 

cos <fi = -RT(a^ai-|-a2a^~),sin <fi = iK{a\a\ — aia\) 1 (6) 

where the constant K = K\Ki- Using the conservation 
character of the square summary, ^sin 2 (f> + cos 2 = 1, 
it is easy to obtain K\K 2 = 1/^/2 (2nin2 + n\ + n 2 )- 
The expectation value and the variance of cos 4> are 

/ a\ (aiai+a2at) , - 

(cos 4>) = =, (7) 

yj2 {2nin 2 + n\ + n 2 ) 
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and 

A(cos <f>) — (cos 2 — (cos^) 2 , 
respectively. Here, the expectation value for cos^ 
_ 1 , <(a+ ai ) 2 ) + <(a 2 4) 2 ) 



( 



cos 



2 (2nin 2 + Tlx + n 2 ) 



(8) 



(9) 



If the matter waves located at two neighboring sites have 
a well-defined phase difference cj>, the variance of cos(/> 
vanishes. Otherwise (e.g., the values of 4> are random), 
the expectation value for cos <f> and the corresponding 
variance will approach zero and 1/2, respectively. The 
random distribution of the phase difference means 'melt- 
ing' of the coherence between the matter-wave clouds lo- 
calized at the neighboring sites. 



2 0.25 




0.25 



FIG. 5: (colour online) Quantum fluctuations for the vortex 
states vs. the ratio U /T for different values of N. 

For the q- vortices (5) with the charge L = 1, by adia- 
batically increasing U/T we calculate numerically the ex- 
pectation values (cos (ft) and variances A(cos0) between 
the sites 1 and 2 (where cos0 = cos(02 — 4>i)) for differ- 
ent U/T starting from the linear limit. For the commen- 
surate cases, the expectation value (cos 4>) grows with 
N and decreases with U/T, but the variance A(cos</>) 
decreases with N and increases with U/T, as shown in 
Fig. 5. These results mean that quantum fluctuations 
decrease with N but increase with U/T. In the limit of 
strong tunneling, i.e. for U/T <C I, the expectation val- 
ues (cos0) — > cos(27r/3) and variances A(cosc/>) — > for 
N — > oo. The limit of strong tunneling with large total 
particle numbers corresponds to the classical limit with 
well-defined phases. The opposite occurs in the limit of 
weak tunneling (U /T 3> I) where the expectation values 
(coscj)) — > and variances A(cos0) — * 1/2 for arbitrary 
N. This means that the distribution of the phase differ- 
ences becomes random, and the classical discrete vortices 
melt completely under the action of strong quantum fluc- 
tuations. Varying U/T < 1 to U/T ^> 1, we recover the 
crossover regime that connects these two limits. 



We suggest that these effects can be studied experi- 
mentally by loading an 87 Rb condensate into the triple- 
well potential formed by a superposition of a harmonic 
potential (with frequency u> ~ 300 x 27rHz) and a Kagome 
lattice [ljj, see Fig- !• The Kagome lattice can be formed 
by laser beams with wavelengths of 1064 nm. The ratio 
T/U can then be adjusted by varying the laser intensity 
about from 1.5 to 15 W/cm 2 to observe the crossover 
regime and the melting of the vortex phase. 

In conclusion, we have studied a quantum analog of 
discrete vortices as the states of interacting bosons with 
a nontrivial phase structure. We have introduced the 
q-vortices for the simplest case of a three-site Bose- 
Hubbard ring and analyzed the effect of quantum fluctu- 
ations on these states. We have found that the melting 
of discrete vortices via quantum fluctuations accompa- 
nies the crossover from superfluid to Mott insulator. We 
believe our findings may initiate experimental efforts to 
observe quantum discrete vortices in Bose-Einstein con- 
densates in optical lattices. 
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